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ABSTRACT 

Broadband optical and IR spectral energy distributions are determined for spectro- 
scopically confirmed z > 2 Lyman break objects in the Hubble Deep Field (HDF). These 
photometric data are compared to spectral synthesis models which take into account 
the effects of metallicity and of internal reddening due to dust. It is found that, on av- 
erage, Lyman break objects are shrouded in enough dust (typically E(B — V) ~ 0.3) to 
suppress their UV fluxes by a factor of more than 10. The dust-corrected star formation 
rate in a typical HDF Lyman break object is ~ 60/i^QQMQyr -1 (qo = 0.5). Furthermore, 
these objects are dominated by very young (<^ 0.2 Gyr, and a median of ~ 25 Myr) stel- 
lar populations, suggesting that star formation at high redshift is episodic rather than 
continuous. Typically, these star formation episodes produce ~ 1O 9 /i^ o 2 o M of stars, or 
~ ^ of the stellar mass of a present-day L* galaxy. 



1. INTRODUCTION 

Recent observational studies of the high-redshift Universe have begun to yield direct informa- 
tion on the nature of the population of "normal" galaxies at z > 2. Color-selected searches have 



started to produce growing, spectroscopically-confirmed samples (e.g., Steidel et al. 1996a; Bteidel 



et al. 1996b| ; Pascarelle et al. 1996| ; Lowenthal et al. 1997; [Steidel et al. 1997| ) , while larger samples 



of objects identified solely on the basis of multicolor photometry have been used to boost statistics 
(e.g., Madau et al. 1996| ; Sawicki, Lin & Yee 1997; Connolly et al. 1997| ). Whether spectroscop- 



ically confirmed or not, such high-z objects are often called "Lyman-break galaxies" or "U-band 
drop-outs". The names reflect the fact that in these galaxies the UV flux is severly attenuated 
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Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. 

2 Based on observations made at the Kitt Peak National Observatory, National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 
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by the combined effects of internal and intergalactic Lyman line and continuum opacity. Unlike 
serendipitously discovered single objects, these new samples of high- 2: galaxies are large enough to 
allow the study of global properties of the z > 2 galaxy population, such as their number density, 
luminosity density, luminosity function, and rates of star formation and metal ejection. 

Despite this recent progress, there is much that still remains unclear. Based on number densi- 
ties, star formation rates, physical sizes, morphologies and masses inferred from spectral line widths, 
Steidal et al. (1996a) and Giavalisco, Steidel, & Macchetto (1996) argue that Lyman break objects 
are the spheroidal components of present-day luminous galaxies, seen in the act of formation. Using 
diagnostics similar to those employed by Steidel et al. and Giavalisco et al., Lowental et al. (1997) 
favor the interpretation that Lyman break objects are either progenitors of present-day low-mass 
spheroidal galaxies or are starbursting building blocks of more massive galaxies of today (see also 
Colley et al. 1997). 

Estimates of star formation rates in Lyman break galaxies are based on their rest-frame UV 
fluxes (e.g., pteidel et al. 1996a ; Madau et al. 1996; Lowenthal et al. 1997| ). However, even moderate 
amounts of dust may significantly suppress the UV flux and, hence, the inferred star formation rate. 
Meurer et al. (1997) have studied the UV spectral indices of HDF Lyman break galaxies and have 
deduced that, if ages of their stellar populations are ^ 0.3 Gyr, reddening due to dust supresses 
the UV flux fairly strongly (on average, a factor of 15 at 1600 A). Using the same spectral indices, 
but assuming older ages (1 Gyr), Calzetti (1997) found the attenuation a factor of 3 lower than 
that arrived at by Meurer et al. (1997). These two different results, obtained from the same data, 
reflect the fact that both dust and aging have similar effects on the UV shape of the SED and that 
one can mimick the effects of the other. This dust-age degeneracy can be lifted by the addition of 
IR (i.e., rest-optical) data, which are sensitive to age but much less so to dust. 

Ages of star-forming objects at z > 2 are not well known. Steidel et al. (1996a) feel that ages 
younger than a few xW 7 yr are unlikely, arguing that such young ages would imply that a large 
number of galaxies have to be bursting simultaneously. However, apart from the gravitationally 
lensed z = 2.72 OB star absorption-line galaxy cB58 ( Yee et al. 1996| ), for which fits to broadband 
optical and IR colors give the age of the dominant stellar population to be ^ 35 Myr flEllingson 
|et al. 19*96 ), no quantitative age measurements exist for Lyman break galaxies. Constraints on 
ages will lift the dust-age degeneracy, in addition to providing clues to the evolutionary history of 
Lyman break objects. 



The Hubble Deep Field (HDF; Williams et al. 1996 ) is a small but intensly studied area of the 



sky, for which both extremely deep multiwavelength imaging (e.g., Williams et al. 1996 ; Dickinson 
|et al. 1997 ; Hogg et al. 1997; Serjeant et al. 1997; Fomalont et al. 1997) and follow-up spectroscopy 
(e.g., Steidel et al. 1996b; Cohen et al. 1996; Lowenthal et al. 1997j ; Zepf et al. 1997] ) have been 
carried out. To date, 17 spectroscopically confirmed z > 2 Lyman break galaxies have been reported 
in the HDF (Steidel et al. 1996a; Lowenthal et al. 1997). 



In this paper we construct the broadband optical and IR (i.e., rest-frame UV and optical) spec- 
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tral energy distributions (SEDs) of the 17 spectroscopically confirmed Hubble Deep Field Lyman 
break galaxies (§ |2|). These observed SEDs are then fitted with dust-corrected spectral synthesis 
models described in § 3.1. These fits show that light from Lyman break galaxies is attenuated by 
significant amounts of dust (typically E(B — V) ~ 0.3; § 3.2) and is dominated by very young 
0.2Gyr) stellar populations (§ |3.3[ ). Furthermore, Lyman break objects have high rates of star 
formation (typically ~ GOhT^MQ-yr -1 for qo 



W 9 h^ M G 



0! 



0.5; § |3.4|) , and quite low stellar masses (typically 
|3.5|) . Since internal reddening in Lyman break galaxies appears to be significant, 



in § LI we calculate its effects on the recently-published estimates of z > 2 star formation and 
metal ejection rates. In § L2 we speculate on the nature of Lyman break galaxies in view of the 
high star formation rates and young ages of stellar populations which we found earlier. Our results 
are summarized in § II. 



2. DATA 



We study the broadband spectral energy distributions for the 17 spectroscopically confirmed 
z > 2 HDF objects reported by Steidel et al. (1996b; 6 objects^], some with revised redshiftsL 1 ]) 
and by Lowenthal et al. (1997; 11 objects). We measured broadband fluxes using both optical 
( Williams et al. 1996j ) and infrared ( |Dickinson et al. 1997 ) data. Specifically, in the optical we used 



the publicly-available Version 2, 4 x 4 binned, HST images of the HDF; in the IR we employed 
the, also publicly available, Version 1 Kitt Peak 4m IRIM observations. The use of the IR data is 
crucial, as at z ~ 3 the optical data alone give coverage only bluewards of rest-2000 A, while the 
IR extends this to rest-5500 A, thereby spanning the age-sensitive Balmer break. Altogether, the 
available photometry spans 3000 - 22000 A (or 750 - 5500 A at z = 3) in 7 bandpasses (C/300, -B45O) 
^606; -^814) J, H, K s ^, although, as will be explained later, we do not include the C/300 and -B450 
data in our fits. 

We smoothed the HST images to match the poorer, ground-based point spread function (PSF) 
of the IR images. We then performed photometry using the PPP faint galaxy photometry package 



3 Steidel et al. report five z > 2 redshifts. Since, morphologically, C4-06 appears to have two distinct components, 
we have treated it as two objects at a common redshift. 

4 The redshift of the Steidel et al. (1996b) object C5-02 has been revised from z = 2.845 to z = 2.008 (Steidel, 
private communication). The reported redshift of C3-02 may also be erroneous: The observed broadband SED of 
C3-02 (Figure [j]) exhibits a break between the 7gi4 and J filters and not, as would have been expected of the Balmer 
break in a z = 2.775 object, redward of J. We note that the absorption lines of C IV 1550, Fe II 1608, and Al II 1670 
could be misidentified as Ly a 1216, Si II 1260, and O I 1303. In this case the redshift of C3-02 would be 2 = 1.95, 
which is consistent with the photometric redshift z p hot = 2.1 of Sawicki et al. (1997) and with the placement of the 
Balmer break between the Isi4 and J filters. While adopting the revised redshift for C5-02, we retain the published 
value for C3-02. We tested the effects of this possible redshift misidentification on our analysis and found them to 
be small. 



5 U 30 o, B450, Veoe, and 7 814 are used to denote the HST filters F300W, F450W, F606W, and F8UW. 
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(Yee 1991). Colors were measured in "color apertures": for each object, a color aperture was chosen 
so as to maximize the S/N. Photometric calibration was based on the zeropoints provided in the 
STScI-HDF and KPNO-IRIM-HDF web pages ( [Ferguson 1996) ; |Dickinson 199^ ), while conversion 
to AB magnitudes (for the IR data) was done using the F u (0 mag) fluxes of Wemstaker (1981). 
The photometry results are presented in Table [I]. The magnitudes listed are "total magnitudes" 
obtained by correcting the object's magnitude within the color aperture by its growth curve in the 
Veso6 image. A typical total magnitude aperture is 3.2". 



3. COMPARISON WITH SPECTRAL SYNTHESIS MODELS 

3.1. The Models 

For the analysis of the stellar populations of Lyman break galaxies, we generated a set of 
model broadband SEDs constructed in the following way: We started with the Bruzual & Chariot 
(1996) multi-metallicity spectral synthesis models; specifically, we used the KL96 theoretical stellar 
spectra^] with the Salpeter 0.1 < 4jr- < 125 IMF. These Bruzual & Chariot models were attenuated 
with Calzetti's (1997; see also Calzetti et al. 1994) empirical reddening recipe for star- forming 
galaxies^] , covering a range of E(B — V) values. The reddened SEDs were appropriately redshifted 
and then convolved with the instrumental transmission curves (in the case of the HST data) or filter 
transmission curves (for the Kitt Peak IR observations) producing a suite of model colors. For each 
object this suite of model colors contains two possible star formation histories (instantaneous burst 
and constant star formation), 221 ages (0-20 Gyr), three different metallicities (O.O2Z0, O.2Z , 
and I.OZq), and a range of reddening {E(B — V) = 0.0 to 0.5 in steps of 0.02). 

Why were these particular ingredients chosen for the models? Initial mass function: The initial 
mass function, at least in low-redshift starbursts, appears to be independent of environment and is 
consistent with the Salpeter IMF with a high (~ 100-M©) upper mass cutoff ( [Stasihska Sz Leitherer 



1996| see also |Massey et al. 1995| ). We have thus chosen to use the Salpeter 0.1 < < 125 IMF 



available in the Bruzual & Chariot (1996) library. Reddening curve: Since, by virtue of their color 
selection criteria ( [Steidel et al. 1996b| ; [Lowenthal et al. 1997| ), the 17 spectroscopically confirmed 
HDF Lyman break objects are star- forming, the use of the Calzetti (1997) reddening curve for 



6 In general, substantial differences still exist between composite SEDs based on theoretical and empirical stellar 
spectra. We have, however, tested our dust and age results with Bruzual & Chariot solar-metallicity empirical models 
and found that the differences were small for our purposes. 

7 Calzetti (1997) gives: 

F o6s (A) = FoWlO- ' 4 ^-^', 

where 

f [(1.86 -0.48/A)/A-0.1]/A+ 1.73 for 0.63^m < A < 1.0/im 

1 2. 656(-2. 156 + 1.509/A - 0.198/A 2 +0.011/A 3 ) +4.88 for 0.12^m < A < 0.63/xm 
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star-forming regions is appropriate. Unlike the commonly-used LMC and SMC reddening curves 
(Fitzpatrick 1986; Bouchet et al. 1985) which are are derived for stars alone, the Calzetti curve 



holds for star-forming regions; it thus automatically includes such effects as back-scattering and 
geometrical distribution of dust and is therefore much more appropriate for correcting the effects 
of dust on the photometry of star- forming galaxies ( |Gordon, Calzetti, fc Witt 1997 ). Using the 



LMC reddening law would produce extinction values 2.5-3 times lower than that which we find 
(§ ^[l]) with the more appropriate Calzetti extinction curve. Metallicity: Damped Ly-a systems at 



z = 2 — 3 have metallicities Z <J O.IZq ( Pettini et al. 1994 ; Lu et al. 1996 ; Pettini et al. 1997a ); thus 



0.2Zq or 0.02Zq are reasonable values to expect for the metallicity of Lyman break galaxies, while 
I.OZq provides a useful check in case the metallicity of Lyman break objects were significantly 
higher than that of damped Ly-a clouds. Ages: At z ~ 3 the Universe is only 1-2 Gyr old and no 
object within it can be older than that; we have, however, allowed for model ages of up to 20 Gyr 
as a consistency check of the fits. Star formation history: In reality, the most likely star- formation 
scenario is that of a burst in which the star formation rate (SFR) declines with time. The two star 
formation histories used (instantaneous burst and constant SFR) thus provide the bracketing cases 
of the likely star formation histories of these objects. Dust contents and ages of stellar populations 
obtained from these two extremes will then bracket the actual values of age and reddening. 

For high-redshift objects, the neutral hydrogen contained in Ly-a clouds along the line of 



sight provides an extra source of opacity in the UV (e.g., |Madau 1995|) . The amount of extinction 
provided by this mechanism is stochastic and depends on the numbers, redshift distribution, and 
column densities of the Ly-a clouds along the line of sight to the target galaxy. The galaxy's 
UV flux is affected blueward of rest-1216 A and is almost completely extinguished below rest- 
912 A. Whereas this effect is of great usefulness in identifying high-z galaxies, its stochastic nature 
introduces uncertainties below ~ 1200 A in the SED models for individual galaxies. These stochastic 
uncertainties, combined with the fact that the dust reddening curve is not well known below 1200 A, 
led us to decide not to fit the two bluest filters (C/300 and -6450)- Thus all the fits presented in this 
paper have been done using the 5 reddest filters (Vqoq, Ism, J, H, and K s ). 

Metallicity (Z = 0.02ZQ, 0.2Zq, and I.OZq) and star formation history (instantaneous burst 
and constant SFR) were fixed, giving six scenarios. In each of these six scenarios, age since the onset 
of star formation and amount of dust were free parameters. The models and the data were compared 
by means of x 2 minimization. As an illustration, Figure p] shows the fits thus obtained with the 
0.2Zq continuous-SFR models. Instantaneous burst models, as well as those for Z = O.O2Z and 
I.OZq, produced fits of similar quality. Note that the two bluest filters, although not used to 
produce the fit, match the best-fitting model SEDs quite well. 

Because reddening increases towards shorter wavelengths, the part of the SED which contains 
information about dust attenuation is the rest-frame UV. Dust attenuation can, however, be mim- 
icked by aging of the stellar population: as the population ages, massive stars die out and the UV 
flux decreases. For example, below rest-3000A a zero-age SED with E(B — V) = 0.3 looks very 
much like a 100 Myr, constant SFR, dust-free one. The availability of IR data allows the lifting 
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of this degeneracy between age and reddening: at z > 2, IR data are probing the rest-optical part 
of the SED and consequently measure the age by means of the size of the Balmer break and, to 
a lesser degree, the slope of the continuum redward of rest-4000A. Metallicity effects could pose 
a problem because of the age- metallicity degeneracy (e.g., Worthey 1994), but we have monitored 



their effects by performing fits with models of different metallicity (0.02Zq, 0.2Z©, and I.OZq). As 
will be seen later, no overwhelming metallicity effects have been found. 

With increasing redshift, the rest wavelengths probed by each filter move more and more 
blueward (by z ~ 3, the H filter is centered at 4000 A) and consequently less and less information 
is available about the age-sensitive rest-optical part of the SED. Thus, past z = 3, only one filter 
(K s ) is redward of the 4000 A break. Because of this loss of age-sensitive information at higher 
redshifts, we have chosen to split the sample at z = 3. We consider fits for the 11 objects at z < 3 
to be much more reliable than those for the six at z > 3. The results and discussion presented 
below are mainly based on the 11 z < 3 galaxies; we will, however, retain the six z > 3 objects for 
completeness and as a consistency check. 



3.2. Dust Content 

For each object the models described above fit simultaneously the dust content and the age 
of the dominant stellar population. The age and E(B — V) parameters of the best-fitting Q.2Z & 
models are plotted in Figure & with error bars corresponding to 90% confidence limits. The Figure 
illustrates that the Lyman break objects are best fitted by models which require non-zero amounts 
of dust. The median values for both the instantaneous burst and constant SFR O.2Z scenarios 
are E(B — V) = 0.28, producing a factor of 16 attenuation at 1600A. This attenuation value is in 
agreement with the result derived from the UV spectral indices by Meurer et al. (1997) (see also 
Calzetti 1997] ) . Note, however, that Meurer et al. had to assume (correctly, as we shall see in § 3.3 ) 



that the z > 2 galaxies are no older than 0.3 Gyr. 

Just as in the O.2Z case, the 0.02Z Q and I.OZq models reqiure significant amounts of dust. 
The results are summarized in Figure || which shows the distribution of E(B — V) values for the 
six different scenarios. Dust is required for the vast majority of galaxies in all of the six scenarios. 
For sub-solar metallicities (0.02^© and 0.2Z & ) typical values are E(B — V) ~ 0.3 and somewhat 
lower (E(B — V) = 0.15 — 0.2) for the unlikely case of solar metallicity. 

Can dust-free models be securely ruled out? The error bars in Figure ^ are 90% confidence 
limits; hence, at 90% confidence, non-zero amounts of dust are unavoidable in all but one or two 
of the 11 z < 3 objects. The presence of dust is further illustrated in Figure ^j, where, for three 
objects, the best possible dust-free models are compared with those in which dust is allowed. More 
qualitatively, the median reduced x 2 are 1-34 and 1.68 for the instantaneous burst ad constant SFR 
0.2Zq fits with dust; the corresponding dust-free values are 5.04 and 10.1. Models with dust give 
much better fits than those without it. Therefore we conclude that dust is present in Lyman break 
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galaxies. A similar conclusion was reached by Ellingson et al. (1996) for the z = 2.72 galaxy cB58, 
which has relatively high-precision IR photometry. 

The high reddening values obtainded here from broadband photometry of Lyman break galaxies 
are in contrast to the negligible extinctions inferred from chemical abundances in damped Lya 
systems (Lu et al. 1996; Pcttini et al. 1997a). However, while quasar sightlines probe essentially 



random regions of damped Lya galaxies, broadband photometry targets regions of star formation. 
Thus, if the star-forming regions have substantially different dust-to-gas ratios than "random" 
parts of high-redshift galaxies, then the discrepancy between the damped Lya spectroscopic and 
our photometric results may simply reflect the differences in the environments that the two methods 
probe. 

Star formation rates have been derived from the fluxes of either Ha or H/3 for the small number 



of high-redshift objects for which these fluxes have been measured (Bechtold et al. 1997; Pettini 



et al. 1997b|) . In these objects, line-based SFR measurements are generally lower than the dust- 



corrected star-formation rates derived from UV fluxes. This apparent discrepancy can, however, 
be caused by a number of effects other than discrepant dust measurements; these effects include 
leakeage of Lyman-limit photons from H II regions, absorption of Lyman-limit photons by dust, 



and a non-standard IMF (see Bechtold et al. 1997 and references therein). These effects may affect 



the reliability of SFR measurements based on emission lines. Ultimately, concurrent measurements 
of Ha/H/5 ratios will provide a robust and independent estimate of the dust content in Lyman 
break galaxies. 

Based on our analysis, the presence of dust in Lyman break galaxies seems unavoidable with 



E{B — V) ~ 0.3 being typical. As will be discussed in § 4T, these amounts of dust are large enough 



to significantly affect estimates of star formation and metal ejection rates at high redshift. 



3.3. Ages of Stellar Populations 

For each galaxy, we shall date the age of the galaxy's dominant stellar population. It must 
be emphashized that the age of the dominant stellar population does not mean the age of the 
galaxy itself, but rather corresponds to the time since the onset of the most recent major episode 
of star formation. Note that the instantaneous burst and constant SFR models are the limiting- 
case scenarios, while a galaxy's star formation history probably lies somewhere between these two 
extremes. Thus, the instantaneous burst model provides the lower bound, while the constant SFR 
model gives the upper bound on the likely age of the galaxy's current episode of star formation. 

Figure 2 presents the parameters of the best-fitting 0.2Z Q models, while Figure 5 shows the 
distribution of galaxy ages for all six combinations of metallicity and star formation history. As 
these Figures illustrate, in the constant SFR scenario the vast majority of z < 3 objects, for 
which spectral coverage is sufficient in the red, are best fit with models for which the ages of 
dominant stellar populations are less than 0.2 Gyr. In the instantaneous burst scenario, the ages 
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of virtually all the objects, including the six with z > 3, are less than 0.1 Gyr. The median ages 
for the z < 3 objects are 36 Myr and 10 Myr in the O.2Z constant SFR and instantaneous burst 
scenarios, respectively. Although objects at z > 3 have insufficient IR coverage to constrain ages 
reliably, their fits are also consistent with very young stellar populations; for the instantaneous 
burst scenarios these ages are as young as for the z < 3 objects. 

How confidently can we rule out older ages of stellar populations? The error bars in Figure || 
are 90% confidence limits. If we consider the constant SFR model (which is more favorable for long 
ages), we see that of the z < 3 objects all but one must have ages < 1 Gyr at the 90% confidence 
level. A further illustration of the necessity of young ages is given in Figure ^: the (forced) 1 Gyr- 
old models produce poorer fits than the unforced younger models. More qualitatively, whereas the 
median reduced x 2 are 1-34 and 1.68 for the 0.2Zq instantaneous burst and constant SFR models 
with free age, the corresponding values for the forced lGyr fits are 355.2 and 6.67, respectively. 
We therefore conclude that the z > 2 Lyman break galaxies are dominated by very young stellar 
populations. 

A possible concern is that the IMF employed to generate the models (i.e., the Salpeter 
0.1 < < 125 IMF) may not reflect the true IMF in the Lyman break galaxies. In particu- 
lar, the determination of the age of the stellar population is primarily sensitive to the size of the 
Balmer break, which grows as the population ages. Since low-mass stars are responsible for the 
Balmer break, an IMF deficient in low-mass stars would produce a Balmer break that is smaller 
than expected for a given age, thereby mimicking the signature of a younger stellar population. 
Consequently, for an IMF deficient in low-mass star, our fits would underestimate the ages of stellar 
populations. The multi-metallicity SEDs of Bruzual & Chariot (1996) are available only for the 
0.1 < -jj- < 125 IMF. Hence, to estimate the impact of an IMF with a high lower-mass cutoff, we 
have compared the predicted colors of Bruzual & Chariot (1993) solar metallicity models with dif- 
ferent IMF mass ranges. Specifically, we have compared the 0.1 < < 125 and 2.5 < -jj- < 125 
Salpeter IMFs. The comparison reveals that, for the young ages in question and if the IMF is 
deficient in low-mass stars, the age of the star-forming population can be underestimated by a 
factor of ~ 5 in the instantaneous burst scenarios; in the constant SFR models the underestimate 
is only at the 10% level. These effects would then increase the median age in either scenario to no 
more than ~ 50 Myr. IMF effects are thus unlikely to affect our conclusion that the Lyman break 
objects have very young stellar populations. 

Figure [7| shows, as a function of redshift, ages of stellar populations obtained with the O.2Z 
SEDs. The age of the qo = 0.5, to = 12-5 Gyr Universe and the age of an object which formed at 
Zf = 4.5 are shown for comparison. Lyman break galaxies are, by and large, dominated by very 
young stellar populations. The absence of older (age > 0.2 Gyr) objects, particularly at the lower 
end of the redshift range covered, is remarkable. Possible reasons for this absence will be discussed 



in 8 4.2. 
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3.4. Star Formation Rates 

To estimate the star formation rates in Lyman break galaxies, we use the constant SFR fits. 
For each of the observed objects, we obtain the SFR by comparing the normalization of the fit to 
that of a lM & yr~ l model of identical age and reddening. For the z < 3 objects, 0.2Z@ constant 
SFR models, and q = 0.5, the median SFR is 59/if o 2 o M yr- 1 {mh^M Q yr- 1 for q = 0.05). The 
distribution of star formation rates in our sample is shown in Figure || 

As a check on the star formation rates derived above from SED fitting, we can use the predicted 
rest-frame 1500 A flux of a 9 Myr old, 1O 6 M galaxy with a 100 < < 1 Salpeter IMF QLeither et 



aL 1995| ). Comparing the dust-corrected, observed fluxes of the HDF Lyman break objects against 
the fiducial flux of Leitherer et al. gives a median SFR of 63/i lo 2 o M0yr~ 1 (go = 0.5), thereby 
confirming the SFR value obtained earlier from SED fitting. 

For illustrative purposes, we have also computed the SFRs under the assumption that the 
Lyman break objects are dust-free. The resulting median star formation rate is 2hiQ M ( r ) yr~ 1 for 
qo = 0.5 (6/ij~ 2 MQyr _1 for qo = 0.05), in agreement with the dust-free SFR estimates given by 
Steidel et al. (1996a) and Lowenthal et al. (1996). However, since models which include dust produce 
much better fits, we prefer the ~ QQh^QQM Q yr~ l (~ 17O/i^ o 2 o M0yr _1 for qo = 0.05) dust-corrected 
rates of star formation. 



3.5. Stellar Masses 

To estimate the stellar masses produced in Lyman break galaxies, we can use the results of 
either the instantaneous burst fits or the constant SFR fits. For the instantaneous burst fits, we 
compute the stellar mass formed in each of the observed objects by comparing the normalization 
of the fit to that of a 1M model of identical age and reddening. For the z < 3 objects, O.2Z 
instantaneous burst models, and qo = 0.5(0.05), the median mass in stars is 1(2) x 10 9 /i^ qMq, or 
about (js) the stellar mass of a present-day L* galaxy (3 x 10 10 /i 10 1 Mq; |Cowie et al. 1995| ). 

Under the constant SFR assumption, the typical duration of a starforing espisode is 2x the 
median observed age. Using the ages and SFRs from the constant SFR fits, we get the median 
stellar mass to be 2(4) x 10 9 /ij~ 2 M o for q = 0.5(0.05) — in agreement with the value derived from 
instantaneous burst fits. As was the case for star formation rates, comparing dust-corrected UV 
fluxes to Leitherer et al. models confirms the above masses. The distributions of stellar masses 
for both the instantaneous burst and constant SFR models are shown in Figure ^(b-c). In both 
scenarios, the stellar masses of Lyman break objects are generally smaller than the stellar mass of 
a present-day L* galaxy. Stellar masses appear to be relatively insensitive to the details of star 
formation history. 

As an illustration, the above analysis can also be applied to the fits obtained with dust-free 
models. For go = 0.5(0.05), under the instantaneous burst assumption, the resulting median stellar 
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mass is 1(4) x 10 9 /iJ" 2 M o . Using the ages and star formation rates of the constant-SFR fits gives 
the median stellar mass of 6(16) x 10 9 /ij~ 2 M Q . These values, obtained from dust-free model fits, 
are surprisingly similar to the stellar masses derived from models which account for dust. In the 
instantaneous burst case, this agreement arises because the dust-free fits give older ages, which 
necessitates larger stellar masses to account for the observed fluxes. In the constant SFR case the 
agreement occurs because the lower star formation rate of the dust-free models are offset by their 
older ages. 

The relative insensitivity of the total stellar mass formed to the assumed star formation history 
and dust content of the galaxy allows us to obtain a robust result for the typical mass produced by 
a star formation episode in a Lyman break galaxy. Averaging the median masses calculated from 
the dust-corrected instantaneous burst and constant SFR fits, we have that the median stellar mass 
produced in an episode of star formation in an HDF Lyman break galaxy is ~ 1.5(3) x 1O 9 /i7/ o 2 o M 
for qo = 0.5(0.05). This median stellar mass is within the 1O 1O M upper limits on the gravitational 
masses of HDF Lyman break galaxies which were obtained by Lowenthal et al. (1997) on the basis 
of Lya emission line widths and is equivalent to ~ ^ (jq) the stellar mass of a present-day L* 
galaxy. 



4. IMPLICATIONS 

4.1. Dust Corrections 

Intervening dust has the effect of absorbing UV photons , thereby suppressing the observed 
UV flux. If, as the SED fits discussed in § |3^ indicate, z > 2 objects have substantial amounts of 
internal reddening, recent measurements of the rates of star formation ([Steidel et al. 1996a| ; |Madau 



et al. 1996; Lowenthal et al. 1997| ) and metal ejection (Madau et al. 1996; fgawicki et al. 199^ ) need 



to be revised. 

Steidel et al. (1996a) have used rest-1500A fluxes to estimate the average star formation rate of 
their Lyman break galaxy sample to be 2/ij~ o 2 o M yr _1 for q$ = 0.5 (and 6/ij~ o 2 o M yr _1 for go = 0.05). 
If we assume that their objects suffer from the same amount of reddening as is typical for the Lyman 
break objects in the Hubble Deep Field (E(B — V) ~ 0.28), then their SFR estimates need to be 
adjusted upward to SSh^MQ-yr^ 1 for q = 0.5 (and to lll/ij~ o 2 o M yr -1 for q = 0.05). 

On the basis of 1500A fluxes, Lowenthal et al. (1997) made estimates of star formation rates 
both in individual z ~ 3 HDF objects and in the volume-averaged population. Reddening correction 
would bring the Lowenthal et al. range of star formation rates to 13 — 35/i{oo^©y r 1 f° r Qo = 0.5 
(31 - lllh^ M Q yr- 1 for q = 0.05). The volume -averaged star formation rates given by Lowenthal 
et al. need likewise be adjusted upwards by a factor of ~ 18. 

Madau et al. (1996) studied a sample of U- and B-band dropouts in the Hubble Deep Field. 
They used the comoving UV (1620 A at (z) = 2.75 and 1630 A at (z) = 4.0) luminosity densities 
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to estimate the volume- averaged rates of metal ejection and star formation. If we apply the typical 
E(B — V) = 0.28 extinction found to occur in z > 2 galaxies, the Madau et al. estimates of both 
the star formation and metal ejection rates need to be revised upward by a factor of 16. This value 
is consistent with the factor of 15 derived by Meurer et al. (1997) on the basis of optical data alone 
and with the (as we have seen, correct) assumption of young ages. Note, however, that the Madau 
et al. (1996) estimates do not include incompleteness corrections, so they remain lower limits even 
after having been corrected for dust. 

As part of their photometric redshift study of the Hubble Deep Field, Sawicki et al. (1997) 
have calculated the expected metal density of the Universe as a function of redshift. To do so, they 
have used rest-3000A luminosity densities. Applying E(B — V) = 0.28 of extinction yields a factor 
of 7 increase in flux at 3000A, requireing a factor of 7 increase in the 2 <^ z ^ 3.5 metallicity of the 
Universe shown in Fig. 10. of Sawicki et al. 

In summary, the values of star formation and metal ejection rates presented in the literature 
are underestimates, since they do not account for intrinsic reddening. Though the exact correction 
factors depend on the wavelength used for the original estimate, once correction for dust is made, 
both star formation and metal ejection rates go up by about an order of magnitude. Note that 
although the total mass of stars formed per object is similar for the dust-free and dust-included 
models (§ [3J^), the difference in the volume-averaged star formation and metal ejection rates of the 
Universe enters via the very different duty cycle of visibility which the Lyman break objects have 
in these two models. 

4.2. Ages and Star Formation Histories of Lyman Break Galaxies 

For all but one of the z < 3 HDF Lyman break galaxies, the ages of the dominant stellar 
populations are < 0.2 Gyr. Since our redshift range (2 < z < 3.5) spans a Gyr in time, the lack of 
a significant number of galaxies with old (> 0.2 Gyr) stellar populations suggests that, in Lyman 
break galaxies, the duration of star formation is short rather than extended and continuous. 

In the rest of this section we will speculate on some of the possible interpretations of this 
episodic star formation. Using simple calculations, we will show that under this single-burst scenario 
the number densities of Lyman break objects, in addition to their stellar masses, are consistent with 
Lyman break objects being progenitors of present-day L < L* galaxies. They are unlikely to be 
the present-day L > L* galaxies caught in the act of creation unless one invokes more complicated 
mechanisms such as mergers or recurrent bursts of star formation. 

In the simplest scenario, akin to the low-mass starburst model of Lowenthal et al. (1997), each 
Lyman break galaxy undergoes a single episode of star formation. Since the stellar mass produced 
in a typical Lyman break galaxy is ~ IO^H^Mq (§ |3.5| ), such a burst would result in a low-mass 
object whose stellar mass is ~ that of a present-day L* galaxy. This typical mass is thus too 
low to account for the whole stellar mass of the bulge of a present-day L* galaxy. 
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The young ages of their stellar populations, together with the fact that they will fade rapidly 
after star formation ends (by ~ 1 magnitude in a mere 30 Myr), mean that only a small fraction 
of such short-burst galaxies are visible at any one time. The redshift range covered by the HDF 
Lyman break sample (2 < z < 3.5), spans a time interval of ~ 0.8 Gyr in a h\oo = 0.7, qo = 0.5 
universe (~ 2.3 Gyr for hioo = 0.7, qo = 0.05). The median age from the 0.2Z & , constant SFR fits 
is ~ 35 Myr. Assuming that a typical star-forming burst lasts twice that time, and that it takes 
~ 30 Myr for the galaxy to fade ~ 1 magnitude (and hence out of the spectroscopic sample), we 
estimate that a typical Lyman break galaxy stays in the sample for ~ 100 Myr. Therefore, for 
every galaxy that is detected, there will be ~ 8 (~ 23 for qo = 0.05) objects which have faded 
beyond spectroscopic detectibility. Hence, in the volume and redshift range sampled, there should 
be a total of ~ 136 (~ 391, for qo = 0.05) visible plus faded objects which are undergoing, or have 
undergone, brief but intense bursts of star formation. 

The comoving number density of this combined (both visible and faded) population is 7 x 
10~ 2 /i 3 00 Mpc -3 (4 x 10~ 2 ft,f 00 Mpc~ 3 for q = 0.05). In the i^-band (which corresponds to ~ rest- 
V), the HDF Lyman break galaxies span a range of ~ 2.5 magnitudes. Locally, over a corresponding 
magnitude range, one expects to see ~ 2 x 10~ 2 ft,^ QMpc~ 3 galaxies with L m (Loveday et 



ah 1992 ; Lin et al. 1996| ). The number density of HDF Lyman break objects is thus similar to, or 
perhaps somewhat higher than, that of local ^qL* galaxies. Thus, on the basis of their number 
densities and stellar masses, Lyman break objects are consistent with being progenitors of present- 
day sub-L* galaxies. 

One could postulate that present-day L > L* galaxies formed through bursts of star formation 
which are analogous to those which seem to be producing the above-mentioned low-mass objects, 
but which are more intense though less numerous than the bursts seen in the HDF. Based on the 
present-day luminosity function of Lin et al. (1996), and assuming that the 17 objects seen in the 
HDF are precursors of present-day sub-L* galaxies, one would expect the HDF to yield ~ 5 Lyman 
break objects with stellar masses > 3 x 1O 1O M0, the expected stellar mass of an L* galaxy. The fact 
that such massive stellar populations are absent (Figure ||) implies that most present-day massive 
galaxies did not form in single bursts of star formation at z > 2. 

There are, however, alternative ways to assemble massive objects. One option is through 
galaxy-galaxy merging: the present-day merger rate is rather low but increases as (1 + z)~ 3-4 , at 



least to moderate redshift (e.g., [Yee k, Ellingson 1995| ; [Patton et al. 1997|) . If this redshift trend 



continues to z ~ 3, then the merging of a dozen or so Lyman break objects can easily result in the 
spheroid of a present-day L* galaxy. 

Another possibility is akin to the "Christmas tree" model of Lowenthal et al. (1997; see also 
Collcy et al. 1996j ), In this model, star formation is episodic but recurrent within each galaxy: star- 



forming episodes are separated by quiescent intervals during which the galaxy temporarily fades 
out of the sample. The underlying stellar population, resulting from previous star-forming bursts, 



would be too faded to be detectable in the presence of an ongoing burst (e.g., Ellingson et al. 1996). 



-13- 



Each of the recurrent star-forming bursts would then add, typically, ~ 10 9 /i 10 qMq of new stars to 
the galaxy, gradually building up its stellar mass. A spheroid containing ~3x W w h^ M G worth 
of stars could be thus accumulated in a dozen or two such recurrent starburst episodes. 



5. CONCLUSIONS 

We have fitted the broadband spectral energy distributions of spectroscopically confirmed 
z > 2 HDF Lyman break objects, using the multi-metallicity spectral synthesis models of Bruzual 
& Chariot (1996). In the fits we have included correction for internal dust extinction typical of 
local star-forming galaxies (Calzetti 1997). The fits also assume that the IMF at high redshift is 
not unlike that seen locally, although some leeway is allowed, particularly in the lower mass cutoff. 

We find that Lyman break galaxies are dominated by very young stellar populations (< 0.2 
Gyr). The absence of objects with old stellar populations, particularly at lower redshifts (z ^ 2), 
implies that star formation in a typical Lyman break galaxy cannot go on continuously for a 
prolonged time. Instead, star formation must occur in bursts of short duration (t < 0.2 Gyr). 

A typical Lyman break galaxy is shrouded in enough dust to suppress its UV flux by a factor of 
~ 16 at 1600A. Consequently, recent UV-based estimates of the rates of star formation and metal 



ejection at high redshift (Steidel et al. 1996a; Madau et al. 1996; Sawicki et al. 1997; Lowenthal et 



al. 1997) need to be adjusted upwards by a similar factor. 

Star formation rates in Lyman break galaxies are high (median of 59/i~ 2 M^yr^ 1 for qo = 0.5). 
These star formation rates typically produce 10 9 /i~ 2 M of stellar mass during the lifetime of a star- 
forming episode; this number is robust as it is relatively insensitive to the details of star formation 
history and dust content. This median stellar mass is equivalent to the stellar mass contained 

in a present-day L* galaxy. Objects with larger stellar masses may be built up through recurrent 
episodes of star formation, or can be assembled through mergers. 

A stellar population will fade rapidly after the end of a brief episode of star formation (~ 
1 magnitude at rest-2000A in the first 30 Myr). Because star formation in Lyman break galaxies 
is of brief duration, a substantial population of objects which have just recently faded out of the 
spectroscopic sample may exist at high redshift. We estimate that if a typical Lyman break galaxy 
stays in our sample for 0.1 Gyr, then there are 8-23 such faded galaxies for every visible one. The 
number density of the combined visible and faded population (a few xl0~ 2 /if o niag -1 Mpc~ 3 ) is 
comparable to the number density of present-day sub-L* {Mb ^ —18) galaxies. 

On the basis of the evidence presented above, we conclude that Lyman break objects are likely 
direct progenitors of present-day sub-L* galaxies, or that they may form luminous galaxies through 
merging or by repeated episodes of star formation. They do not appear to be steadily star-forming 
direct progenitors of present-day massive galaxies. 
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Table 1. Photometry 
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a Since object C4-06 has two very distinct components, we have treated it as two objects at the same redshift. C4-06a 
is the lower object in the C4-06 panel of Fig. 2 of Stcidcl ct al. (1996b), while C4-06b is the upper object. 

Note. — Upper limits (2a) are indicated with a ">" . Errors are indicated in parentheses: they are la uncertainties 
expressed as percentages of the flux. 

References. — (S) Stcidcl ct al. (1996b); (L) Lowcnthal ct al. (1997) 
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Fig. 1. — The best-fitting O.2Z , constant SFR models. Dust content and age are free parameters. 
Although their fluxes are plotted, the t/300 and -B450 data were not used in the fit (see text). The 
model SEDs shown include the high-z Lyman suppression of the UV flux. The bottom two rows 
contain objects at z > 3 — objects for which fit quality suffers from the fact that only one filter is 
present redward the Balmer break. Instantaneous burst and O.O2Z and l.OZ models produced 
fits of similar quality as the ones shown here. 



-20- 



0.4 



0.2 







0.4 



0.2 



- 



constant SFR, O.2Z 



error barfl are 90% 
confidence limits 



I I I M il 



instantaneous burst, O.2Z 



error barfl are 90% 
confidence limits 



0.01 



0.1 

age 



10 



Fig. 2. — Reddening and age of the best-fit 
Z = 0.2Z & models. The top panel shows the 
constant SFR fits and the bottom one is for the 
instantaneous burst model. Age is the time since 
the onset of star formation. Galaxies at z < 3 
are shown as solid symbols while those at z > 3, 
for which fit quality is poorer, use broken ones. 
Error bars come from % 2 fitting and correspond 
to 90% confidence limits. Z = 0.02Z Q and l.OZ 
models produce similar results. 
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Fig. 3. — Intrinsic reddening for the 6 models 
of different metallicity and star formation his- 
tory. The shaded histogram corresponds to those 
galaxies at z < 3, while the unshaded one is for 
all 17 objects (i.e., including the lower-quality 
z > 3 fits). Arrows indicate the median values 
of E(B — V) for the z < 3 galaxies. Substantial 
amounts of dust are present under all scenarios. 
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Fig. 4. — Examples of fits with and without dust 
for the Z = 0.2Z@ constant SFR models. The 
broken line shows the best-fit models for which 
age was allowed to vary freely, but E(B — V) was 
held at 0. The solid line shows the best-fit model 
(as in Figure |l|) for which both age and dust 
were free parameters. Models with dust produce 
better fits than the dust-free ones. 
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Fig. 5. — Ages since the onset of the cur- 
rent episodes of star formation. The shaded his- 
togram corresponds to those galaxies at z < 3, 
while the unshaded one is for all 17 objects. Ages 
> 0.5 Gyr are grouped together in the right-most 
bins (separated from the other bins by the dot- 
ted lines). Most Lyman break galaxies are dom- 
inated by very young stellar populations. 
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Fig. 6. — Examples of fits with and without 
forced age for the Z = 0.2Z & constant SFR mod- 
els. The broken lines are the fits for which age 
was fixed at 1 Gyr (while reddening was allowed 
to vary freely). The solid lines are the fits where, 
as in Figure [l], both age and reddening were 
free parameters. Instantaneous burst models of 
1 Gyr age would produce an even stronger dis- 
agreement with the data, since in those models 
the UV-producing massive stars would not have 
continued to be replenished as in the constant 
SFR models. 
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Fig. 7. — The ages of dominant stellar popula- 
tions of Lyman break galaxies determined by fit- 
ting 0.2Zq model SEDs. Error bars correspond 
to 90% confidence limits. For comparison, the 
dotted lines show the age of the universe and 
the age of an object which formed at zj = 4.5; 
both are for a fi = 1 universe with a present- 
day age fixed at to = 12.5 Gyr; note, however, 
that the stellar population ages do not depend 
on the value of Hq. Objects at z > 3 suffer 
from poor coverage above ~ 4000 A and, con- 
sequently, have age estimates of lower quality. 
Models with Z = O.O2Z and I.OZq produce 
ages similar to these shown here. Stellar popu- 
lations of the majority of z > 2 HDF galaxies 
appear to have undergone recent (t < 0.2 Gyr) 
episodes of star formation. 
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Fig. 8. — Distribution of star formation rates 
obtained from the 0.2Z Q constant SFR fits. A 
qo = 0.5 universe was assumed. The shaded his- 
togram is for the eleven z < 3 galaxies, and the 
unshaded one includes all 17 objects. The me- 
dian value for the z < 3 subset is indicated. 
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Fig. 9. — Stellar masses from the 0.2Z© fits. 
Shaded histograms are for z < 3 galaxies, and 
the unshaded ones include all 17 objects. Arrows 
indicate median values for the z < 3 subset. A 
qo = 0.5 universe was assumed. Panel (a) shows 
the distribution of stellar masses obtained from 
the instantaneous burst fits. The stellar mass of 
a present-day L* galaxy (3 x 1O 1o /i5~ o 1 o M ; Cowie 
et al., 1995) is indicated for comparison. Panel 
(b) shows the distribution of stellar masses de- 
rived by multiplying the objects' star formation 
rates and the ages of their stellar populations. 
Under both the instantaneous burst and con- 
stant SFR scenarios, stellar masses produced by 
the observed episodes of star formation are gen- 
erally substantially less than the stellar mass of 
a present-day L* galaxy. 



